What are the interactions in quantum glasses? 
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The form of the low-temperature interactions between defects in neutral glasses is reconsidered. 
We analyse the case where the defects can be modelled either as simple 2-level tunneling systems, 
or tunneling rotational impurities. The coupling to strain fields is determined up to 2nd order in 
the displacement field. It is shown that the linear coupling generates not only the usual Ising- 
like interaction between the rotational tunneling defect modes, which cause them to freeze around 
a temperature Tq, but also a random field term. At lower temperatures the inversion symmetric 
tunneling modes are still active - however the coupling of these to the frozen rotational modes, now 
via the 2nd-order coupling to phonons, generates another random field term acting on the inversion 
symmetric modes (as well as shorter-range l/r"' interactions between them). Detailed expressions 
for all these couplings are given. 
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I. INTRODUCTION 

One of the most subtle and peculiar problems in con- 
densed matter physics concerns the nature of the " glass" 
state, and of the associated glass transition. This prob- 
lem (described by P. W. Anderson[lj in 1995 as "the 
deepest and most interesting problem in solid-state the- 
ory") concerns the overwhelmingly dominant compo- 
nent of the physical world as we experience it, ic., non- 
conducting solids that are not ordered in regular crys- 
talline arrays. In fact the glass problem actually involves 
two separate features. One is the remarkable universality 
displayed in the low-T quantum properties 0,3, and the 
other is the high-T behaviour shown in the vicinity of the 
glass transition itself H. 

At first glance it seems implausible that these two fea- 
tures of glasses could be related - they occur at very 
different energy scales. Elsewhere we argue that there 
may be an interesting connection, which depends on cer- 
tain novel features of the interactions in these systems. 
The purpose of the present paper is to investigate the 
form of these interactions in some detail. We derive a 
number of new interaction terms, which arc presented in 
the form of two new effective Hamiltonians for neutral 
glasses, one valid for higher temperatures, the other in 
the low-T limit. 

We begin with a brief introductory review of the 
physics of neutral glasses, particularly in the low T quan- 
tum regime. We then derive the form of the defect- 
phonon interaction terms (section |Tl|, including both a 
direct linear coupling to the lattice displacement field, 
and a coupling to the gradient of this field. We then 
calculate, in sections |lll] and HVl the effective coupling 
between defects induced by these interactions. It is 
shown in section IIIII that the linear coupling not only 
gives the well-known Ising coupling between the rota- 
tional tunneling modes, but also a random field acting 
on these modes. However there is no such linear cou- 
pling between phonons and the 'inversion symmetric de- 
fects' (ones which are symmetric with respect to inver- 
sion about the local lattice site). In section |IV] we cal- 



culate the coupling of these defects to gradients of the 
phonon field, and show that this coupling produces an- 
other weaker coupling to a random field generated by 
the now frozen rotational modes, as well as a short-range 
coupling between the inversion symmetric tunneling de- 
fects. 

The main point of the present paper is to find the 
correct quantitative form of the effective Hamiltonian for 
these systems, after integrating out the phonons. We end 
up with an effective Hamiltonian at high T which involves 
only the rotational defects, and then another quite differ- 
ent low-T Hamiltonian which describes a set of tunneling 
inversion symmetric defects, coupled to each other and 
to the random field generated by the frozen rotational 
modes; this includes a number of new terms. It turns 
out that these results have very interesting implications 
for the physics of glasses, which are discussed in detail in 
Ref.[|. 



A. Universalities in the low-T quantum state 

The first glass conundrum, strongly emphasized by 
LeggcttQ, is the apparent universality in the low-T prop- 
erties of a huge variety of disordered systems below a 
temperature Tq ~ 1 — 3 K, regardless of the amount of 
disorder. The most striking universalities are seen in 

(i) the Q-factor Q{uj, T) for torsional oscillations of the 
system, which is related to the phonon mean free path 
l{ijj,T) and the phonon thermal wavelength A(a;,T) by 
Q = 27r//A. One finds Q that Q{uj,T) shows a pro- 
nounced plateau for T < Tq (down to a lower tem- 
perature which decreases with uj), and that the value 
Q = Qo ^ 600 varies over only a factor ~ 2 — 3 between 
many different materials, even though the intrinsic disor- 
der (measured by e.g., defect concentration x) may vary 
over several orders of magnitude. 

(ii) The " Berret-Mcissner ratio" c^/c|| between longi- 
tudinal and transverse sound velocities in the same low- 
T regime. A remarkable linear relationship is foundQ 
between c± and C|| , for a variety of materials includ- 
ing amorphous oxide, semiconducting, polymer, metallic. 
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and electrolyte glasses, in which c\\ varies by a factor of 
5. 

These are only the most striking of the low T uni- 
versalities - there are others In recent years the 
experimental groups of Osheroff [3 |3l , HunklingerQ, and 
EnssfToj have pushed experiments to very low tempera- 
tures [~ 0(1 mK)], and found a host of interesting new 
results, including intrinsic dipolar "hole-digging" in the 
many-body density of states [7|, and a remarkable spin 
coherence phenomenon!^, [lo| which comes from nuclear 
quadrupolar interactions in neutral glasses. The dynami- 
cal hole-digging persists to the lowest temperatures, giv- 
ing ever sharper features in the density of states; it is 
associated with non-exponential relaxation and 'aging' 
behaviour of the dielectric constant of the system. It 
would be of interest to continue these measurements well 
into the regime, if possible. 

Although we are cle arly dealing here with a resonant 
tunneling phenomenon [ll|. q ualitatively similar to that 
in tunneling spin svstems[l2l. [isj the glass problem ap- 
parently involves cooperative tunneling of at least cou- 
pled pairs of tunneling systems [1, [ll|, and the resulting 
low-T "universal" state apparently involves some funda- 
mental new physics. Although a number of theoretical 
scenarios have been pro posed to describe this univer- 
sal physics H, [m, Q m, most of which argue that it 
must involve strong coupling between the relevant low- 
T modes, there is no complete consensus at the present 
time. We note in passing that although the universalities 
occur in the same temperature range as the well-known 
regularities [3| in thermal and transport properties in 
glasses (such as the specific heat Cy(T) ~ AT, or the 
thermal conductivity K{T) ^ BT^), these latter can all 
be understood in terms of the well-known picture of 
non-interacting two-level systems(TLSs). On the other 
hand the dynamics of dipolar hole-digging certainly re- 
quires interactions for its explanation, between what- 
ever modes are exhibiting low-T quantum fluctuations, 
whether these be pairs of TLSs[ll| or some more compli- 
cated set of modes 0, [13]. 



B. The high-T glass transition 



aging effects, in the vicinity of the glass transition (as 
noted above, these effects are also found in the low-T 
quantum regime, both in neutral glasses where univer- 
sality is seenQ and in electronic glasses [isj). 

At first glance there seems to be no obvious relation 
between this high-T behaviour, which is characterised by 
thermally-activated processes of great complexity, and 
the low-T behaviour. A number of attempts have been 
made in the last couple of years to give a general theory of 
the glass transitionQlli,!!^]. Two of them[li,[23| made 
no connection to the low-T regime, instead concentrat- 
ing on the vast number of thermally-activated processes 
coming into play near Tq. The Moore- Yeo theory also 
makes a very interesting connection between the critical 
behaviour of supercooled liquids near Tq and Ising spin 
glasses in a magnetic field. 

However, there are several arguments that indicate 
that there may be a connection between the physics be- 
low Tq and that near Tq- The first two are experimental. 
It was already noticed by Berrct and MeissnerQ that 
the low-T phonon relaxation time Tmin shows systematic 
connection to the values of Tq across the whole range of 
glasses, with Tmm on Tq^. We have already remarked on 
the rough proportionality Tq/Tq ~ 40; and indeed one 
can argue that the coupling 7 between phonons and the 
low-T tunneling entities (whatever they may be) is intrin- 
sically related to both Tq and the phonon velocities [14 1. 

These observations suggest that there may be some 
kind of unified theoretical framework which could de- 
scribe both the high and low-T properties of glasses. 
Such a theory would not only be of major interest (an- 
swering the question posed by Anderson [l|) but it would 
also clearly give us a new blueprint for theories of other 
complex systems. Such a framework has actually been 
proposed very recently by Lubchenko and Wolvnes . 
In this theory the basic objects are "tunneling centers" 
which comprise 200 atomic units, and which can be 
used to describe both the low-T dynamics and the dy- 
namics near Tq- We note that these "tunneling entities" 
are very different from the TLSs that have been used to 
describe many of the low-T experiments (sl. [ol. [Tol. [Til. [l6t . 
although it is argued that their behavior will be quite 
similar [ill . 



In most glasses there is actually a glass transition at 
a temperature Tq much higher than the crossover tem- 
perature Tq to the universal quantum regime (typically 
Tg/Tq ~ 40). There are several universal features of this 
transition as amongst which one may single out 

(i) the characteristic range of relaxation times r in the 
system, and their characteristic T dependence, summa- 
rized in the Vogel-Fulcher scaling law, which shows that 
the value of Tq we use depends on the timescale in inter- 
est; (ii) the "entropy crisis", in the range Tk <T < Tq, 
were Tk is the Kauzman temperature, and where one 
finds a supercooled glass entropy lower than that of the 
crystalline solid; and (iii) the existence of highly non- 
exponential relaxation, and characteristic memory and 



C. Nature of interactions in glasses 

We now come to the question to be addressed by this 
paper. One might think, in view of the generality of 
the phenomena discussed above, that they ought to be 
independent of the detailed nature of the interactions be- 
tween tunneling entities in the low-T regime. Actually it 
is widely assumed in the glass literature that it is enough 
to include only dipolar strain-mediated interaction, with 
the addition of electric dipolar interactions where neces- 
sary. This assumption rests on a number of microscopic 
calculations done over the years, both for interacting 
TLSs[li[, [m [H [H, [13, and for the more comphcated 
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systems of discrete rotators that exists in orientational 
glasses [25! [26I [27l. [28j. These calculations (particularly 
those for orientational glasses) are very lengthy, and for 
this reason have not been attempted with any generality 
except by a few authors. Nevertheless, the conclusion 
has been that the effective interaction between 2 tunnel- 
ing entities a distance R12 apart is 

H,nt - gia^/Ruf (1) 

where g is a coupling ~ 0(1 eV), with both dipolar strain 
and electric dipole contributions. The length ao is a typ- 
ical lattice distance; and we have suppressed an angular 
factor which has roughly dipolar symmetry. We note that 
all of the principal scenarios for the low-T behaviour of 
glasses assume ([1]) to be true; moreover they rely upon 
it in an essential way. 

But is ([1]) really correct? In this paper we shall ar- 
gue that in fact ([T]) is incomplete, and that the correct 
form contains extra terms of some importance. These in- 
clude another term falling off like 1 / i?f 2 j which however 
leads to a random field acting on each tunneling rotator. 
There are also terms which are weaker and which fall off 
faster, which would be less important except that they 
act on 2-level systems that do not see the first random 
field. The net result of this is that we derive two effective 
Hamiltonians for the system, one which is valid at higher 
temperatures around Tq and the other at much lower T, 
apparently around the temperature Tq which defines the 
crossover to the universal properties. 

II. DEFECT-PHONON INTERACTIONS 

In what follows we will be dealing with neutral glasses, 
i.e., we ignore metallic and superconducting glasses. This 
of course still includes the overwhelming majority of ma- 
terials on earth, from rocks and minerals to a galaxy of 
insulating compounds (based mostly on transition met- 
als), along with a huge variety of natural and artificial 
organic systems (including polymers). 

In spite of this variety, and regardless of whether one is 
dealing with a strongly disordered amorphous systems or 
very weakly disordered systems like substitutional elec- 
trolyte glasses, the two interactions of main interest are 
those involving strain fields and phonons, and those in- 
volving the interaction of electric fields with local charge 
distributions. 



A. Simple model for Defects 

Our tactic in this paper will be to start with a toy 
model which describes a class of very simple systems, and 
then argue that the important features of this model can 
be generalised to a much wider variety of glassy systems. 

Consider the situation depicted in Fig. [1] in which we 




(a) (b) 

FIG. 1: (a) Off-center impurity with 4 equivalent sites and (b) 
orientational impurity in the dumbbell approximation, with 2 
equivalent states. To first order in the displacement, the two 
opposite sites in (a) are also equivalent, leading to an effective 
2-state impurity. 



reduce the system to a 2-d plane, and look at defects in 
this plane. The underlying local symmetry of the sys- 
tem is assumed to be of square plaquettes, and the de- 
fects can either be substitutional defects, able to occupy 
one of 4 states in the plaquette, or orientational defects, 
able to rotate between 4 orientational states. Under cer- 
tain circumstances, to be discussed below, we can make 
the 'dumbbell approximation, in which defects rotated 
by 180° are considered to be indistinguishable - we then 
treat these 2 states as identical, and the 4-state system 
reduces to a 2-state system. 

Now if the concentration of these defects is low we can 
assume that they do not disturb the underlying lattice 
symmetry, and interactions between 2 defects, even if 
they are distantly separated, will be between 2 plaquettes 
which are oriented along the same axes. More generally, 
when the defect concentration is much higher, one may 
have the situation shown in Fig. [21 where the 2 plaquettes 
are slightly distorted, and also rotated with respect to 
each other. 

It might be objected that the situation depicted in Fig. 
[2] is not very realistic, in that a high defect concentra- 
tion would so severely disrupt the square symmetry that 
the plaquettes themselves would not only be rotated, but 
that their shapes would be severely distorted, so that 
no clear local lattice structure could be defined. Actu- 
ally this is not the case - quite surprisingly, the situation 
in even rather strongly amorphous glasses does not con- 
form to the common caricature in which they look like 
frozen liquids. Instead, at short length scales the un- 
derlying lattice structure is still quite recognizable, and 
the system more resembles a 'frozen liquid' of very small 
polvcrvstals[29| (actually the instantaneous state of quite 
a few liquids also looks like this!). 

A good tutorial example of systems like our toy model 
is provided by substitutional electrolyte glasses, where 
the defects can be characterized very precisely. Canonical 
examples are KCli_xLix, KCli_x(OH)x, KBri_x(CN)x, 
and so on[23, HO]. Many experiments in these systems 
have been done in the very dilute regime, with concen- 
trations in the range 10~^ < x < 5 x 10^^, where one 
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FIG. 2: Two nearby impurities, with distorted plaquettes. 
Filled small circles denoted occupied impurity sites. 



can ignore the interactions between the widely spaced 
impurities [Sll. [s^. [ssj . However there are also many stud- 
ies of concentrations up to x ~ 0.1 or even higher. When 
X > 2 X 10~^, interactions between the impurities clcaijy 
dominate the physics in most of these systems[30, 33, 
[35j . This is shown in, e.g., the saturation in the dielectric 
response for higher concentrations (ssj . 

In real 3-dimensional systems like (KCl)i_xLix, 
(KCl)i„x(OH)x, or (KBR)i^x(CN)x, one has either off 
center point-like impurity states (as in (KCl)i_xLix, 
with 8 available positions for the Li impurity inside a 
given KCl lattice 'cage'); or else orientational (as in 
(KBR)i^x(CN)x, where the CN impurity can lie in one 
of the 8 directions along (1,1,1) and its equivalents). 
Previous treatments of the microscopic interactions in 
these systems have taken one of two routes. Michel and 
collaborators [2^123 have used a detailed microscopic de- 
scription of the orientational and translational degrees of 
freedom (and the coupling between the two of these) with 
the goal of describing systems where orientational impu- 
rities like CN can rotate inside their host 'cages'. Sethna 
and collaboratorslll,!!!,!!^! have used a somewhat more 
general phenomenological description in terms of TLS 
impurities. 

The great advantage of beginning with these simple 
systems is that one can set up a well-controlled theory 
beginning with the case of dilute defect concentration x. 
There are good arguments, based on the experiments on 
electrolyte glasses 13^, 'sTUsaJslJsl [H and other amor- 
phous systems ^3. Z. 8, .oTIlH' 1 1 . 16l|. that the theory must, 
when X is larger, flow towards the strong coupling regime 
(and perhaps to some universal quantum regime) . But it 
makes no sense trying to explore the strong-coupling limit 
until the form of the interactions has been established in 
weak coupling. The key question of interest here (viz., 
what are these interactions?) can only rigorously be ad- 
dressed by starting from a system where x^ 1. Later 
we will argue that our main results will survive well into 
the universal regime, and so are of much more general 
applicability. 

Returning now to Fig. [l] we divide the defects into 
two groups, viz. (i) those in which there is an inver- 
sion symmetry relating pairs of states, where inversion 
is made with respect to the relevant lattice position; and 



(ii) those where there is no such symmetry. In the square 
plaquette system we see that both simple impurities like 
Li and orientational impurities like CN or OH can be 
in one 4 different states, related by 90° rotations of the 
plaquette, and there is no inversion symmetry. However 
in some cases, one can treat the states related by 180° 
as physically equivalent (either because they really are 
equivalent, or because at the energy scale of interest the 
difference is unimportant), and then we can assume in- 
version symmetry. In this case the state space on the pla- 
quette is 2-dimensional, with states oriented along one or 
other of the lattice diagonals. 

To describe the system, the simplest representation is 
just a 4-state one in which the system can hop between 
any of the 4 plaquette sites. However we will use a slightly 
different one, where we begin by defining a set of oper- 
ators {fi} such that rf flips the state on the plaquette 
through 180°, and another set of operators {Si} such that 
Sf rotates the state on the plaquette through 90° in ei- 
ther direction. This representation is used because we 
will see that the impurity-lattice interaction turns out 
to depend crucially on whether inversion symmetry is 
obeyed, and so this way of setting up the description al- 
lows us to distinguish between operations which are or 
are not invariant under inversion symmetry. 

In what follows we will first be considering the effect of 
phonons on the rotational tunneling degrees of freedom. 
In general, before we take interactions with phonons into 
account, these will be described by a Hamiltonian 



H. 



def 



(2) 



where Dj is a tunneling matrix element and Ej is any 
stray field acting on the rotational defect variable Sj . At 
lower T, once the rotational degrees of freedom are frozen 
out, only the inversion tunneling processes are left, and 
these are described by a bare Hamiltonian 



H. 



(r) 



dof 



A -f^ 



(3) 



where again an "external" effective longitudinal field 
acting on the j-th dipole is allowed. 

Often, in the study of the possible phases of the system, 
the high-T tunneling is dropped since it is typically too 
small to influence the nature of the phase, nor its stability 
to effective random fields [11]. We will keep it simply 
because we wish to trace how the effective Hamiltonian 
evolves as we lower the energy scale. 



B. Coupled Defect-Phonon System 

Consider now the coupling to the phonon modes in the 
system. The Hamiltonian for this system is given by: 



H ^ H, 



def 



ph 



iiit 



(4) 
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where fldcf is the bare defect term just discussed, and 
the phonon system is described by 

Here P and X represent momentum and displacement 
operators for phonon modes of wave vector q and branch 
^, and M is the mass of the elementary cell of the 
medium. 

We split the defect-phonon interaction iJint into two 
terms, writing 

ffint = ^^i"* + H'f (6) 

where the lowest order term is linear in the phonon dis- 
placement field: 

a. 13 j ■' / 

(7) 

and the next term is a non-linear coupling of the defects 
to the gradient of this field, of form 

""'-^p-M'.^^rJ^ (S) 

a, 13,5 J jt J 

The linear coupling term H™^ contains first a 'volume 
coupling' with coefficient 77, which is independent of the 
defect position or orientation, arising because the defect 
has a different volume from the host and this locally 
strains it. Then there is the usual interaction 7Q/3(rj) 
between the defect and the strain field in the "TLS" or 
"dumbbell" approximation, which changes sign with S^, 
since the phonon displacement fields are sensitive to the 
defect orientation. The size of these interactions is not 
easy to calculate - however 7^^ (r j ) can be measured, and 
is on a typical scale 7 ~ 1 eV. Estimates of 77 give similar 
numbers but the ratio 77/7 must certainly vary from one 
system to another (calculations arc difficult because the 
volume change and charge redistributions caused by the 
defect will interact with each other). 

We notice that when inversion symmetry is preserved, 
there is no linear interaction between the fj and the 
phonon field. The 2nd term H™'^ arises when we relax 
the assumption of inversion symmetry, so states rotated 
by 180° arc distinguishable (either because we relax the 
"dumbbell" approximation, and consider the dipole mo- 
ment of the impurity or molecule, or because we con- 
sider tunneling point impurities). In this case the effec- 
tive 'dipole' represented by the difference between these 
states can interact with the gradient of the phonon dis- 
placement field, in the form given in H™^ . However this 
interaction is much weaker, and has rarely been consid- 
ered before. Certainly we do not expect it to affect the 
physics near the glass freezing temperature. However we 
will see that this interaction is important at lower tem- 
peratures, where other degrees of freedom freeze out. At 



the present time we cannot give more than a rough esti- 
mate for the size of the coefficient C,"'^^ of this interaction; 
noting that its dimensions are 

p[^a/35] = EL (ie., energy 
times length) , and that the characteristic energy scale of 
defect interaction energies \s Uo ^ ^ eV, and the char- 
acteristic length scale ad of defect dynamics is roughly 
1 1, we would expect that |C"''''| ^ I eV x X. To get 
an energy from this we need to divide by the typical dis- 
tances between impurities, given by Id ~ aoX~^/^ , where 
Oo is a lattice length. Thus we expect that the energy 
scale associated with this coupling is less than 7 and ?/ 
by a factor adX^I'^ j ao- From now on wc will assume that 
this energy scale associated with Q is considerably smaller 
than 7, 77. 



III. INTER-DEFECT EFFECTIVE 
HAMILTONIAN: HIGH ENERGIES 

Let us begin by assuming that we are at a sufficiently 
high energy scale that we can neglect the weaker interac- 
tion ; we are then concerned with a set of phonons in- 
teracting with the \Sj\ TLS variables, ignoring the {tj} 
variables. We wish, in this approximation, to calculate 
the effective Hamiltonian of the system, initially just up 
to 2nd-order in the defect-phonon coupling. This is an 
old problem, but we shall see that even here there are 
new things to be discovered. Let's look first at the simple 
problem of 2 interacting defects in this dumbbell approx- 
imation. We split the interaction term into two, to 
understand the effect of the orientation and volume parts 
separately, and consider 2 defects at positions ri and V2- 
Then the system has bare Hamiltonian: 

If we now integrate out the phonons we will generate, at 
lowest order in terms of Ising form (proportional to 
I1I2SIS2), cross terms giving a local field, proportional 
to 777(5*^ -I- 51), plus an energy shift ~ 771772. The Ising 
form has been known for a long time. The second term, 
when summed over all spins apart from one given spin, 
simply leads to a random field at the site of that spin - 
this term is not usually considered. In this section we 
first sketch the derivation of the Ising term, primarily 
to establish notation, and then derive the random field 
term. 
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A. Ising interaction term 

We define the Fourier transformation to momentum 
space as 



9,M 



TV 



(10) 



where Sq^^^a is a phonon polarisation index. Then we 
have 

1^7 = --^I^Tla/J^^q.Meg.^.aJg/je''"^? + (1^2). 

(11) 

To find the interaction in 2nd-order perturbation the- 
ory we minimize the potential energy, i.e. the sum of the 
second term in Eq. ([5]) plus the interaction term. Straight- 
forward calculation results in energy terms proportional 
to {Sf)"^, {^i)"^ and the Ising interaction term of inter- 
est, proportional to S'lS'f. Let us define the notation 
R-12 = xi — X2, and use an acoustic approximation, in 
which the longitudinal phonon frequency cOqi = ciq, the 
transverse phonon frequency is U!q± = c±q. We also have 
the identities: 







SM=ii,-L2 eqM"eq^/3 - q^qp/q^- (12) 

Then the Ising interaction = '^fl 'S'i'S'2! where 
^12 = - 2^ MA4, ,2 1^ eq^ceq^iqpqsli li 



X COS [g(a;i - ^2)]. (13) 
Summing over polarization indices gives 



U- 



12 



E 



1 \ ^ Qi3 7(5 



NM 



1 1 



qgqiiq^qs 

q' 



COS (q • R12) 



-1 a/3 aS 

—f 2^ 2^ • ^12) (14) 



NM 



a, 0,5 -L 

which when the sum over momenta is performed, gives 
the real space form 



TTZZ _ 9l2 



(15) 



with the interaction 



5?2(n) 



~1 „a/3„.a5 



^ 2" ) 7r''72''(- [SapSjS + ^ajSps + SasSpj] 



+ 3 [SpsnaTi.^ + S-ysnanp] — Ibuanpn-^ns) (16) 



Here the reduced variable n = R/|R|, and we suppress 
the "{12}" subscript on R12 and ni2 to keep things un- 
cluttered. This interaction has a rather complicated an- 
gular dependence, coming both from the anisotropy of 
the medium, and from the that of the strain interaction 
lap (r) . If we assume a completely isotropic medium with 
degenerate longitudinal and transverse phonons, and also 
make the simplification of anisotropic 70/3 (r), so that 
7a/3(r) —* ^oSazSpz, we get the strictly dipolar form 



[3cos2 6l(n) - 1] 



(17) 



where 6'(n) is the angle between the unit radius vector n 
and the z-axis; the characteristic coupling energy go = 
il"^ / pc^) is now evident. 

As noted before, these results are well-known, and were 
first derived [2ll [22| in the 1970's. The essential result is 
that one has derived an effective Ising interaction; since 
the sites of the defects are random, the tunneling terms 
make this system behave as a quantum Ising model with 
random interactions, and Hamiltonian 



(18) 



where the tunneling amplitudes {Di} are typically much 
smaller than the nearest-neighbour Ising interactions. 
However it turns out that this Ising interaction is not 
the only term that is important. 



B. Random Field term 

We now include the cross-terms oc 777, coming from 
2nd-order perturbation theory in the interaction H™^ in 
Eq.®. Using similar manouevres as in the calculation 
above, one now obtains another term in the effective in- 
teraction of form = ^12 ("^1 + '^'2)7 where 

C^12 (R12) = -T7T7 >^ — ^ > — ^ COS (q • R12). 



a/3 



q 



(19) 

We call this a random field term because if we take a 
given spin in the system, say 5*^, and then sum the inter- 
action Ul^ (Ry ) between Si and the volume terms coming 
from all the other defect sites {_?'}, we get a field h\ acting 
on SI at site i which varies from site to site in a random 
way, because of the random positions and orientations of 
the arguments R^ . Note, incidentally, that this effective 
random field interaction has contribution only from the 
longitudinal phonons. 

This interaction is important since, as we show else- 
where, it actually destroys the bulk glass transition^. 
Let us now explicitly derive its form. To do this we first 
evaluate the tensor 



l-ali 



E 



9a 9/3 , n^ 

COS (g • R). 



(20) 
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Changing the sum to an integral we get: 
V 



I a 13 



{2-KRf 



a q — 5— cos iq ■ n) 

q 



V 



{2nRy 



(21) 

where again n = R/|R|. From symmetry we have I^p = 
o-5ai3 + ^ntjii^. Consider first X^Q^aa- ^''^^ hand 

this sum can be written as / d'^gcos (9 • n) = 0. On the 
other hand the sum equals 3a + h, leading to the identity 
3a + 6 = 0. 

Similarly, the scalar 

^/„^n„n^ = j d'qM^I^ cosiq-n). (22) 

In calculating the right side we can choose n in any direc- 
tion, say in the z direction. This leads to the expression 



27r / dqrq,- / dq^ ^ ^ ^ cos q^ 



and then using 



= 1 - 



(23) 



(24) 



and straightforward integration we find that the integral 
equals — 47r^. on the other hand, X)a/3 IctfS^a^fS = a + b, 
leading to the identity a + b = — 47r^, and together with 
the above identity (3a + 6 = 0) to the result that a ~ 
27r^, b = — 67r^, and 



lafs ~ — 27r^(3nQn^ — Sap)- 
Thus, one finally obtains 



U 



zO 



(25) 



(26) 



Q/3 



This result shows a much less complicated angular de- 
pendence than the Ising interaction (|16p : in general we 
see that they depend differently on angle, simply because 
the Ising g^^{r) is essentially a dipole-dipole interaction 
whereas this mixed term is a dipole-monopole interac- 
tion. 

If we make the isotropic assumption that "fap{r) — > 70, 
then UI2 0. On the other hand if 7^^ = j^z, then 
we get a dipolar interaction: 



g^°(n) 

i?3 



where 



5^°(n) 



[3 cos^6l(n) 



(27) 



(28) 



where now the interaction energy is (72^77/471^0^). In 
any case, both this mixed term and the Ising term end 
up having the same 1 / R^j spatial form, but their charac- 
teristic energies are different. 



Thus, as a result of the added volume term in the 
impurity-lattice interaction, the effective Hamiltonian in 
the " dumbbell" approximation is not the quantum Ising 
model (jl8p . but the quantum random field Ising model, 
with Hamiltonian 



Hrf{{Sj}) 



\ JJZZ QZ QZ 



B,Sl) (29) 



where Bi = J^j ^if ■ This random field has mean zero, 
as by symmetry there is no preferred direction. Its typ- 
ical value is given by Bp ~ rj^yi/ {AiTpcf), where x is 
the concentration, since it is dictated by the typical dis- 
tance between nearest impurities. Thus, the typical size 
of both interactions is rather similar; we expect that 
^ I S j t^i^ I ('y/l) • typical glassy systems this 
means that they are both ~ xJ7o, where fJo ~ 1 eV. This 
means that unless x < 10~^, the typical size of these ran- 
dom fields \Bi\ 3> I?o, where Do is a typical tunneling 
amplitude. Note, that the standard deviation of the dis- 
tribution of the random fields, B ijj^/x/ (47rpc^) Bq, 
as it is dominated by the rare events of pairs of impurities 
occupying nearest neighbor lattice points (sec Ref. [ioj . 
noting the trivial relations between high and low impu- 
rity densities). 

The Random field Quantum Ising Hamiltonian (j29p 
has very different properties from the simple quantum 
Ising system - apart from anything else, the random field 
can actually destroy the glass transition^. We do not go 
into these questions here, but note instead that the most 
important effect of the Ising interaction term, for all but 
very dilute glasses, is to freeze the tunneling of the Sj 
variables, except for a very small fraction ~ Do/xUo of 
systems that happen to be on resonance. 

However this is not the end of the story at all. This is 
because, as noted above, the tj variables do not have a 
linear coupling to the phonons, and so to linear order they 
experience neither impurity-impurity interactions nor a 
random field, and thus they are still free variables. This 
is why we now have to go to the higher coupling terms. 



IV. INTER-DEFECT INTERACTIONS AT LOW 
ENERGY 

Let us now go to an energy scale very much less than 
the putative glass freezing temperature. Now we start 
from a Hamiltonian given by 



E7, 

3,01. (3 



a, 13,5 j 



Q/3i5 9^Xja 



(30) 



where the tunneling amplitude describes the 180° fiip 
of defects. The first term in this effective Hamiltonian 
is just the random field Quantum Ising model derived in 
(pg)) above. 
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Now let us write an approximate version of ([50]) , which 
takes account of the fact that the set of spins {S^} have 
almost entirely frozen into some random configuration, 
with expectation values {(Sj)}, because of the strong 
Ising interaction between them (but without long-range 
glassy order) . Now in this approximation we can simply 
ignore the dynamics of the {Sj} entirely, and replace ((30|) 
with another Hamiltonian, which is approximately valid 
when T <^ Tq, given by: 



3 a, 13,5 j 



This Hamiltonian is only valid to the extent that we can 
ignore any tunneling of resonant {Sj} variables. 

Consider now the interaction terms in Eq. (j3ip . If we 
now integrate over the phonons, one expects the interac- 
tion between fj and the phonons [the 3rd term in (j30p ] to 
give an interaction ~ Q^jT^Tj. But we notice that there 
will also be two cross-terms between the two couplings to 
the phonons, giving an interaction ~ liCj^i (Sj) and an 
interaction ~ ■ When summed over the sites {j}, 

both of these terms must give random fields acting on the 
Tj variable. The first of these random fields comes from 
the frozen Sj degrees of freedom, behaving as a quenched 
impurity distribution coupling to the degrees of free- 
dom. The second just comes from summing over all the 
scalar volume distortions from these frozen impurities. 



which come into the calculation of this interaction when 
we have only 2 impurities. Using Eg. ljlOp we obtain 



2 ■ 



(33) 



Again, minimizing the sum of Vst and the potential term 
in Eq. ([5]) one obtains the interaction energy between the 
two impurities, given by: 



(31) U^2^-Y1 



qti 



(34) 

Notice that this interaction is odd in R! This is be- 
cause the first derivative term is imaginary, and the sec- 
ond derivative term is real, leading to a sin (q • R) inter- 
action. Physically, this is clear since an impurity variable 
which is odd under inversion symmetry, like r2, must have 
an interaction which is odd in R with other impurities, 
if these are either substitutional (ie., causing scalar per- 
turbations), or are even under inversion. 

Using again the acoustic approximation, i.e. uoqi = 
ciq,u!qj_ ~ c±q and the identities in (|12p . we have an 
interaction 



TfSr _ _ 



1 



(35) 



A. Interactions involving the {fj} variables 



This is a complicated integral, because of the large num- 
ber of different components of momentum involved. To 
evaluate it we begin by writing it in the form: 



Since 7, 77 3> Ci it follows that the two random 
field terms acting on the {tj} interaction will be much 
stronger than the Ising interactions between them (quite 
different from what occurs for the {Sj}). We therefore 
deal with the two random field terms first, and then look 
at the Ising interaction between the {tJ }. 



1. 



The SiTj interaction 



tjSt _ 
^12 — 



1 



NM 



where we have define the 5th and 3rd rank tensors 

qgqpq'fqsqr, sin (g ■ R^ 

Q 



(36) 



(37) 



The S^Tj interaction is not only typically the largest 
(although 7 and rj are of the same order, 7 is usually a 
little larger), but also the most tedious to calculate. To 
do this we begin by considering the 2 interaction terms 



Vi 



St 
2 



EQ/3 



8x1/3 



27 



^ ^ dx2sdx2ri'^' 



7,(5,7; 



(32) 



Sl3Sn - E 



qpqsqr, sin (q • R) 



q^ 



(38) 



Let us start by calculating Gpsij- As a first step we 
write it in the form 



V 1 



G 



(39) 
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where we have defined the integral 



where 



^isSti = I a q — 



(40) 



Already here we see that the interaction has a spatial 
dependence ~ In fact, any interaction involving 

the Tz impurity will have a 1 /R"^ or larger power spatial 
dependence (we find that the r^r^ impurity- impurity in- 
teraction has a dependence), and must therefore be 
treated as a short range interaction in 3D. 

By symmetry the function GfjSri can be written in the 
form 

Gpsrj = A{Spsn,j + 63r,ns + Ss,jnp) + Bupngn,,. (41) 

We wish to find the coefficients A and B. We do this 
along the same lines as in our calculation for the random 
field term of the impurities. First we note that 



P /'^3„ gog/3g7g^g')Sin(gTi) 

i'aP'iSTi ~ J 1 ^4 • 

From symmetry, one can write: 

Fafi-ySn = aF^p^gjj + bF^^p^g^ + cF^p^g^ (50) 

where 

+ (77<->a + 77^/3 + 77^7 + 77<-^(5) 
Fap-ySrj ^ '^Q.^n^naii,, + all distinguishable permutations 



(51) 



We now need to evaluate the three coefficients a, b, and 
c. To do so let us first consider the term 



J2 Gaa„n,, =J2[ '^^llv^v sin (g • n) (42) J2 ^""/S/^')"') = J ^^lil ' ^ (q-n). (52) 



is a scalar. Thus the integral can be taken for t] = z, and 
is zero. On the other hand the sum on the left side is 
5A + B, and therefore we get 



5A + B = 0. 



(43) 



af3rj 

The equality results because the sum gives q'^q^j in the 
numerator. This expression is a scalar, which therefore 
equals J (Pqqz sing^ = 0. A careful evaluation of the sum 
on the left side gives 35a + 146 + c and therefore we obtain 
as a first relation for a,b,c that 



We now consider 



Er< /'j3 qnqsq^^npnsnjj . 
Gi3s^nf3nsnn = 2^ / " g — sm [q ■ n). 

l3Sri pSrj ^ 

(44) 

The left side equals 3A + B. On the other hand, since 
the expression is a scalar the integral can be calculated 
for n = Uz, i.e. 



3 9|sin(g^ 

g „2 



/•OO /"OO / 

2n / dqrqr / dqzqzSui{qz) [l - 

Jo J~oo \ 



ql + ql 



= -127r^ 



(45) 



Thus we find that "iA + B = -12t:'^ , and since 5^H-S = 
we find that A = Gtt^, B — — 307r^, and therefore finally 
we have 

G/35,, = 7r^[6(5/35n,, + 5f),^r\s + Ss^^np) - 30n/3n5n,,] (46) 
and hence 

V 

Gpsjj = 47ri^4 ['^('^/3^"'> + SpjjUs + SsrjUp) - 15n^nA-n,,]. 

(47) 



Let us now calculate Tap^Si-f- We first define 
T ~ ^ ^ F 

■TaPySi-i — ^27r)3 _fj4 "/37i5r; 



(48) 



35a -I- 146 + c = 0. (53) 
To obtain a second relation we look at the scalar: 

,3 qjqsqr,n^nsn,j . 

q— —:r sm [q ■ n) (54) 

g^ 



7(5 r; 



Since the expression is a scalar, we can choose n ~ Hz, 
and the integral becomes J (P qq^ sin qz / q^ — — 127r'^, as 
was calculated above. Summation of the left side gives 
21a + 125 + c, and we thus have a second relation between 
a, b, and c, viz.: 



21a+126 + c = -127r^. 
To get a third relation we consider 



(55) 



^ i^a/375»,nan0n^n5n,, = J d^qql sin qz/q'^, (56) 

where we have used, as above, the fact that the expres- 
sion is a scalar and took n = n^. The integral can be 
calculated using again the identity ql/q'^ = 1 — q^/q^, 
and is found to equal zero. The sum on the left side then 
gives the relation 



12a 105 H- c = 



(57) 
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If we now take these three relations together, we find the 
desired resuhs for a, b and c: 



a = 247rV5 
b= -1387rV5 

-2 , 



10927r75. 



(58) 



which can then be inserted into Eq. (|5D)) to get a final 
result for the interaction in the form given in (j36p . 

Now we can write the form of the effective random field 
term that this leads to, after noting again that almost all 
the {Sj} are frozen, and sum over all sites apart from 
a given site i. Then we must get a term in the low-T 
effective Hamiltonian of form 



(7), 



(59) 



where the random field b\^^ is given by summing over 

'e just derived. Or 



sites in the interaction we have just derived. One gets 

^,(7) ^ ^ 

2tt^NM ^ R ^ ^ 

J ''•I ap-yoT] 



2tt^ 

c 



2 "a7 '-'pSrj 



,(60) 



where the tensor function F^J^^^^^ is just 



3 9Qg/3'?7959r, sin(q- riij) 



= [aKp^s., (n^, ) + bFlu^,^ {n,, ) + cF^^^sr, (n„- )] (61) 

with the 3 tensors in the 2nd form given by substituting 
the unit riy = i?ij/|Rjj| for n in ([5T|) : and the tensor 

^/35?) given by the form in ([46]) after the same substi- 
tution has been made, ic.. 



(62) 

This is the first random field term acting on the {fj}. 
We see that the order of magnitude of the interaction is 
given by 



1^: 



(7)1 



V — 

PC' ^ i?^ 



(63) 



where 7 and C ^-rc typical values of the corresponding 
tensors. 



between the volume interaction at site ri and f|. Simi- 
lar manouevres to the ones used above then lead, in the 
acoustic approximation, to an interaction term 



which reduces to 
-1 



E ^ E ^1^2 ^''^agag,, sin (q • R)t| 



-1 



NM ^ q 



9 '-'a-7 



(65) 



NMcf 



jSrj z 



(66) 

where the definitions of Q and G are the same as above 
[cf. Eqs. (|46)l and (|47)l ]. Note that since the volume term 
couples only to the longitudinal phonons, C/,,t- depends 
only on q. 

The random field resulting from this term is then given 
in the form of an interaction 



where the random field bi is given by 



with G'^g^ given by 



3pc2 ^ ^7^')^^^ 



(67) 



(68) 



above, and where 



E r4 • 



PC- Y 



(69) 



Now we can finally write the random fields acting on 
the {tj}, coming from the gradient phonon term in the 
defect-phonon interaction, in the form of an interaction 
term 



E^^ 



(70) 



where the random field bi is just given by the sum of the 
two contribution we have found, ie.. 



l(7) 



(71) 



where the two contributions are given by eqtns. (j36p and 
([68)) respectively, and where \b^i^ lb^i^\ ^ 7/77. 



2. The volume interaction term 

Now let us consider the other source of random fields 
acting on the {fj}, coming from their interaction with 
the volume term of other impurities. This is calculated 
by considering the cross-term arising from the interaction 



Vr^r 



27 



a 



K2j7 



(64) 



3. Ising interaction between the {fj} 

In the same way we may evaluate the Ising inter- 
action between the {f}. The result is an interaction 
Jij^i^j ^ where the interaction coefficient 



pel R% 



(72) 
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Now we see that not only is the interaction coefficient 
here much smaller than that for the random fields (since 
C ^ '777), but this interaction also falls off much faster 
with Rij, going like l/i?fj. For this reason we do not 
calculate the exact coefficient Jij here: this calculation 
is rather lengthy, and it simply multiplies the right hand 
of ([72|l by a complicated angular factor ^ 0{1). Note 
that even though is so small, it will still affect the 
dynamics of the system at low T, since there will be still 
resonant {fj } that can tunnel. 



B. Effective Hamiltonian: low-Energy Form 

Let us now summarize what we have. After integrating 
out the phonons, we can now say that we have ended up 
with an effective Hamiltonian which, if wc still treat all 
the variables as operators, takes the form 



H, 



ff 



-Y}D,S^ + A.r/] + V^ff (73) 



where the interaction now contains the following terms: 

^e// X! ^if'^i^j + X! ^^^^ + X! ^f/^i'^j 

ij i ij 

+ T.^^^ + llJ^^^^^h (74) 

i ij 

where we have written the terms in decreasing order 
of their strength. The 2 random fields in this effective 
Hamiltonian arise from the phonon-mediated coupling of 
the {5?} and the {ff } variables to the volume distortion 
caused by the defects. 

The above effective Hamiltonian must be used if we 
want to analyse the dynamics of these variables. At high 
energies ~ Tq, we can entirely drop all the terms involv- 
ing the {ff} in the interaction Ve//, since these inter- 
actions are all ~ 0(C) and too weak to play a role. If 
we ignore the small number of spins Sj that are in reso- 
nance, then we can also assume that the {S^} variables 
are frozen by the strong Ising interaction U^^' . As dis- 
cussed elsewhere, the effect of the random field Bi is then 
to destroy long-range glassy order @. 

Now suppose we go to low energy scales. If wc continue 
to ignore the small quantum fluctuations in the expec- 
tation values {(iS'j)} brought about by tunneling of the 
small concentration of resonant Sj , then we can treat this 
distribution as frozen. Then the low-energy Hamiltonian 
simplifies very considerably. We get 



H, 



eff 



(75) 



where now the interaction term has the much simpler 
form 



eff 



in which the random field bi was calculated in the last 
section, and, as noted before, it is much larger than the 
Ising interaction Jij. 



V. SUMMARY AND REMARKS 

The purpose of the present paper was to give a detailed 
treatment of the phonon-mediated interactions which ex- 
ist in a neutral glass, taking into account not only the 
usual linear coupling between defects and the phonon 
displacement field, but also the coupling to the gradient 
of the phonon field. The net result of this was that wc 
found an effective Hamiltonian for the two sets of vari- 
ables (the rotational tunneling variables {Sj} and the 
inversion tunneling variables {fj}), which contained tun- 
neling terms for each, Ising interaction terms for each, 
and various effective random fields which act on both 
variables. We note that similar issues arise in quantum 
spin glasses [3^ lU where random field terms also 
have a profound effect (just as they do in classical spin 
glasses [43); however real spins are in many ways quite 
different from defects, and this means that there is no 
simple relation between the two systems 

These calculations were all done in the framework of 
2nd-order perturbation theory in the interactions, and 
they were all done assuming that the background lat- 
tice could still be meaningfully defined, at least in local 
"patches" around each defect. Thus at first glance the 
calculations here are only rigourously valid for low defect 
concentration x. Two possible problems then arise at 
higher x. First, one might object that the 'patch' picture 
must eventually break down - we have argued in section 
|TT]that this is not the case, because even rather strongly 
disordered glasses still do have local crystalline order. 

The second more serious problem is that one expects 
higher-order interactions to come in at higher x and these 
will mix the various interactions wc have derived here. 
This problem of higher-order corrections is notoriously 
difficult, since a hierarchy of lo gs i s generated once one 



39j . We do not attempt 



(76) 



integrates over multiple sites 111, 
to discuss it here, but simply note that our results in- 
evitably changes the results of these higher-order calcu- 
lations, because of the new terms we have found in this 
paper (in fact almost all calculations of these higher or- 
der terms include only the tunneling terms and the Ising 
interactions, without any random fields). Thus we expect 
the results here to have important consequences for the 
discussion of the nature of glasses, and we have developed 
some of these elsewhere @. 
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